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Using a sample of 225.3 million J/tp events collected with the BESIII detector at the BEPCII 
e + e~ collider in 2009, searches for the decays of 77 and 77' —¥ % + e~i> e + c.c. in J/ip —t 4>V and <jrq' 
are performed. The cp signals, which are reconstructed in K + K~ final states, are used to tag 77 
and r\' semileptonic decays. No signals are observed for either r\ or 77', and upper limits at the 90% 
confidence level are determined to be 7.3 x 10~ 4 and 5.0 x 10~ 4 for the ratios B ^Z~7* 7V 1 "e+c- c -) 

and B( - T g^7_ > °4-^- 4 ^' c '' ) , respectively. These are the first upper limit values determined for 77 and 77' 
semileptonic weak decays. 

PACS numbers: 13.20.Gd, 14.40.Be, 13.20.Jf, 12.60.Cn 



I. INTRODUCTION 



Weak decays of quarkonium states such as 77, 77', 
J/ip and T, etc., offer a window into what may lie 
beyond the standard model (SM) [H-[i]. The reason 
for the expected sensitivity is that the rates of the 
quarkonium weak decays are expected to be tiny in the 
framework of the SM Q. As originally pointed out 
by Singer the weak decays 77 — > ir + l~Di (I = e, 
fi, and charge conjugate state implicitly included) are 
purely second class with a vector-type coupling in the 



SM (see Ref. for the definition of the second class 
current), and hence vanish in the limit of exact isospin 
symmetry. They occur in the SM in first order in the 
weak interaction, but only due to G-parity breaking 
effects, i.e. due to electromagnetic corrections and the 
mass-difference of the u- and d-quarks [Icl - [l7l |. For 77 
semileptonic weak decays, a one-loop calculation was 
performed in chiral perturbation theory within the SM, 
including a systematic treatment of the electromagnetic 
contributions to 0(e 2 p 2 ) (e and p are electromagnetic 
coupling and typical momentum transfer in the decay as 
defined in Ref. [l2|), and a rather accurate upper bound 
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for the branching fraction of 77 — > n + l~vi is predicted to 
be 2.6 x 10~ 13 . Therefore, any observation of 77 — > nlvi 
violating this bound would be a clear indication for new 
physics beyond the SM. 

The decays 77 — > nlvi can be used to probe some types 
of possible new charged current interactions [l(J [Hi • A 
rather old suggestion would be the introduction of a new 
second class vector current for the r\ — > it transition [l5| . 
Scalar-type charged current four-fermion interactions can 
arise in gauge theories for example from the exchange 
of charged Higgs bosons in the two-Higgs-doublct 
model [IH HH. Also light leptoquarks [20j |. occurring 
naturally in grand unified theories and composite 
models, may enhance the 77 — > ixlvi branching faction 
considerably [2lj|. For example, by considering scalar 
or vector type interaction, the branching fraction of 
77 —} n + l~vi was estimated to be 1CT 8 - 10" 9 [H |23|, 
which is a few order of magnitudes higher than that 
in the SM. Therefore, searches for the 77 — > n + l~Di 
and rj' — > t: + 1~Di at the branching fractions level of 
1CP 8 — 1CP 9 and below will provide information on the 
new physics beyond the SM. At present there is no 
experimental information on the decays rj —> ixlvi. In this 
paper, we present measurements of branching fractions 
of 77 and rj 1 — > n + e~D e decays. This analysis is based on 
(225.3±2.8) x 10 6 J/ip events [H[, accumulated with the 
Beijing Spectrometer III (BESIII) detector [Hj], at the 
Beijing Electron Positron Collider II (BEPCII). 



II. THE BESIII EXPERIMENT AND MC 
SIMULATION 

BEPCII/BESIII HI is a major upgrade of the BESII 
experiment at the BEPC accelerator. The design peak 
luminosity of the double-ring e + e~ collider, BEPCII, 
is 10 33 cm" 2 s _1 at the center-of-mass energy of 3770 
McV. The BESIII detector has a geometrical acceptance 
of 93% of 4-7T and consists of four main components: 
(1) a small-celled, helium-based main drift chamber 
(MDC) with 43 layers, which provides measurements of 
ionization energy loss (dE/dx). The average single wire 
resolution is 135 fim, and the momentum resolution for 
charged particles with momenta of 1 GeV/c in a 1 T 
magnetic field is 0.5%; (2) an electromagnetic calorimeter 
(EMC) made of 6240 Csl (Tl) crystals arranged in a 
cylindrical shape (barrel) plus two end-caps. For 1.0 GeV 
photons, the energy resolution is 2.5% in the barrel and 
5% in the end-caps, and the position resolution is 6 mm in 
the barrel and 9 mm in the end-caps; (3) a time-of-flight 
system (TOF) for particle identification (PID) composed 
of a barrel part made of two layers with 88 pieces of 5 cm 
thick and 2.4 m long plastic scintillators in each layer, 
and two end-caps with 96 fan-shaped, 5 cm thick, plastic 
scintillators in each end-cap. The time resolution is 80 ps 
in the barrel, and 110 ps in the end-caps, corresponding 
to a 2cr K/tt separation for momenta up to about 1.0 
GeV/c; (4) a muon chamber system made of 1000 m 2 



of resistive-plate-chambers arranged in 9 layers in the 
barrel and 8 layers in the end-caps and incorporated in 
the return iron of the super-conducting magnet. The 
position resolution is about 2 cm. 

The optimization of the event selection and the 
estimation of physics backgrounds are performed using 
Monte Carlo (MC) simulated data sam ples . The 
GEANT4- based simulation software BOOST [26j | includes 
the geometric and material description of the BESIII 
detectors, the detector response and digitization models, 
as well as the track records of the detector running 
conditions and performance. The production of the 
J/t/j resonance is simulated by the MC event generator 
kkmc [271; the known decay modes are generated by 
EVTGEN [28| with branching ratios taken from the 
Particle Data Group (PDG) tables [2!| and determined 
by the Lundcharm model lundcharm (3(| for the 
remaining unknown decays. The analysis is performed in 
the framework of the BESIII offline software system [3lJ 
which takes care of the detector calibration, event 
reconstruction and data storage. 



III. DATA ANALYSIS 

A. Analyses for ij and rj' — > -K + e~v e 

In order to detect rj and rj' — > n + e~i? e decays, we 
use J/ip — > (f>rj and <prj' decays. These two-body decays 
provide a very simple event topology, in which the <p 
signals can be reconstructed easily and cleanly decaying 
into K + K~ . The reconstructed 4> particles can be used 
to tag 77 and 77' in order to allow a search for their 
semileptonic decays. In addition, the 77 and 77' decays 
are easy to define in the lab system due to the strong 
boost of the (f> from J/ip decay. 

Charged tracks in the BESIII detector are 
reconstructed using track-induced signals in the 
MDC. We select tracks within ±10 cm of the interaction 
point in the beam direction and within 1 cm in the 
plane perpendicular to the beam direction. The tracks 
must be within the MDC fiducial volume, | cos(9| < 0.93 
(9 is the polar angle with respect to the e + beam 
direction). Candidate events require four charged tracks 
with net charge zero. The TOF and dE/dx information 
are combined to form PID confidence levels for the 
7r, K, and e hypotheses; each track is assigned to the 
particle type that corresponds to the hypothesis with the 
highest confidence level. To suppress background from 
J/tp — > 4>rj (rj'), where 77 (77') decays into nonleptonic 
modes, the electron candidate is further identified with 
the ratio of deposited energy in the EMC to track 
momentum, E/p, which must be larger than 0.8. We 
further require that E/p should be less than 0.8 for the 
pion candidate to suppress background from J/tp — > <prj 
(rj "/e + e~) decay. 

Showers identified as photon candidates must satisfy 
fiducial and shower-quality requirements. The minimum 
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energy is 25 MeV for EMC barrel showers (| cos#| < 0.8) 
and 50 MeV for end-cap showers (0.86 < | cos6»| < 0.92). 
To eliminate showers produced by charged particles, a 
photon must be isolated from any charged track by more 
than 20° if not specified otherwise. 

Since the mass of neutrino is almost zero and it is 
invisible in the detectors, a one-constraint (1C) kinematic 
fit is performed to constrain the missing mass of the 
reconstructed tracks to be zero, and \ic < 200 is 
required. The 1C fit improves the resolution of recoil 
mass of the K + K~ system by a factor of 2.5 for 77 case 
or a factor of 1.6 for 77' case. After the 1C fit, the 
missing momentum P m ; ss = |P m i ss | can be calculated; 

here, P m iss = — {P<f> + Ptt+ + P e ~) m the rest frame of 
J/ip, and we require that the missing momentum should 
be larger than 0.03 GcV/c to suppress backgrounds from 
final states with only four tracks, such as J/ip —> <Ptt + tt~ 
(<p -> K + K~). We count the number iV showor of EMC 
showers that could originate from a Kl or a photon, 
and require that iV s hower be zero in the region inside a 
cone of 0.3 (1.5) rad around the direction of the missing 
momentum for J/ip — > <pr](r)') [r](j]') — > 7r + e~t> e ]. These 
requirements reject most 77 and rj decays into nonlcptonic 
final states. They also eliminate most backgrounds from 
multi-body decays of J/ip — > </> + anything. The different 
requirements on the cone angle for the r\ and 77' cases 
are made because of the following two reasons: firstly, in 
the J/ip — > <f>r](r]') decays, the booster for 77 is stranger 
than that for 77' in the central of mass energy of J/ip, 
which leads to a larger open angle for the 77' decay 
products than that for the 77 decay products in the 
detector. Secondly, the most dangerous backgrounds are 
from 77(77') — > 7r + 7r~7 decay, in which one of the charged 
pions is mis-identified as an electron. Meanwhile, the 
decay rate for 77' — > 7r + 7r~7 is more than 6 times larger 
than the rate for 77 — > tt + tt~j [29j . 

Figures [1] (a) and (b) show the invariant mass 
distribution of K + K~ candidates, tuk+k--: after the 
above selections. Clear <fi signals are seen. The 
invariant mass of ir + e~V e can be obtained as 777 7r + e - Pe = 

{E n+ + E e - + E p ) 2 - (P n+ + P e - + P p ) 2 , where E„ = 

. = iPnissI and Py = P miss . Figures [2] (a) and (b) 
show the 777 w + e - Pe distributions for events with 1.01 < 
m K+K - < 1.03 GeV/c 2 for the decays J/ip — > (p>r\ (77 — > 
ir + e~v e ) and J/ip — > <pr( (77' — > n + e~v e ), respectively. 
No events are observed in the 77 and 77' signal regions. 
The signal regions for 77 and 77' are defined in the ranges 
[0.51, 0.58] and [0.92, 0.99] GcV/c 2 , respectively, on the 
mass of TT + e~D e . 

We use MC simulated events to determine selection 
efficiencies for the signal channels and study possible 
backgrounds. With phase space MC simulations, we 
obtain efficiencies of 17.9% and 17.4% for 77 and 77' 
scmilcptonic decays, respectively. According to the study 
of the J/ip inclusive MC sample, more than 20 exclusive 
decay modes are identified as potential background 
modes, and are studied with full MC simulations in 
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order to understand the backgrounds. The sources of 
backgrounds are divided into two classes. Class I: The 
background is from J/ip — > <f>rj(rj'), <p K + K~ and 77 
(77') decays into other modes than the signal final states. 
We find that the expected number of background events 
from this class is 0.18 ± 0.05 (0.58 ± 0.09) in the signal 
region for 77 (77'). Class II: It is mainly from J/tp decays to 
the final states without 77 or 77', such as <fiir + ir~ , 0/ o (98O) 
(/o(980) -)■ 7T+7T-), and K*°K*° (K*° -» K^). The 
expected number of background events from class II is 
0.05 ± 0.04 (0.45 ± 0.13) in the signal region for the 77 
(77') case. The total number of background events is 
0.23 ± 0.06 (1.03 ± 0.16) in the signal region for 77 (7/). 
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FIG. 1: The m K + K - distributions of candidate events: (a) 
for J/ip — > (pVi (b) for J/tp — > <pvj ' . The arrows on the plots 
indicate the signal region of cp candidates. 
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FIG. 2: The m^+e-p^ distributions of candidate events: 
(a) for J/ip — > (p-q (77 — ► Tx + e~u e ); (b) for J/tp — > (prf 
(r/ —5- -K + e~i> e ). For both (a) and (b): the data (dots 
with error bars) are compared to the signal MC samples 
(red dashed histogram) and the expected backgrounds (solid 
blue histogram). The arrows on the plots indicate the signal 
regions of 77 and 7/ candidates. 



After all selection criteria are applied, no event 
survives in the 77 and 77' signal regions. The signal 
components and the expected background shapes are 
projected and compared to data for both 77 and 77' cases, 
as shown in Figs. [2] (a) and (b). We set an upper limit 
at the 90% confidence level (C.L.) to be N* h = 2.36 

(A^' L = 1.59) for 77(77'), using the POLE++ program [H[ 
with the Feldman- Cousins frequcntist approach [33j. The 
information used to obtain the upper limit includes the 
number of observed events in the signal region, and 
the expected number of background events and their 
uncertainty. 
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Analyses for 77 (77') — > tv + tv 7r°(r/) 



The 77 — >• ir + ir~ir° and 77' — > ir + Tr~r] decays are 
also studied in J/tp — > <prj and <firj' processes, in order 
to obtain the ratio of 6(77(77') — ► n + e~v e + c.c.) to 
6(77 — > 7r + 7T _ 7r°) (B(r) — > tt + tt^7])). The advantage of 
measuring the ratios of semileptonic weak decays over 

B{r)' — nr^ e~ y e +c.c.) 



hadronic decays B ^^L r ^-X^)' ) and 
is that the uncertainties due to the total number of J/tp 
events, tracking efficiency, PID for kaon and one pion, 
the number of the charged tracks, and residual noise in 
the EMC cancel. 

The selection criteria for the charged tracks are the 
same as those for the J/tp — > 0?? (?/), 77 (77') — > n + e~D e 
decays except for the electron identification requirement. 
The candidate events are required to have two charged 
kaons and two charged pions with opposite charge. 
In addition, two photon candidates are required to 
reconstruct ir a — > 77 and 77 — > 77 in the 77 — > 7r + 7r~7r° and 
7/' — ► 7r + 7T _ 77 decays, respectively. The photon candidates 
are required to be isolated from all charged tracks by 
more than 10° which is different from the selection 
criteria for the J/tp — > (pV (77'), 77 (77') — > 7r + e~;/ e decays in 
order to improve the efficiency of TT°(r]) reconstruction. 
A four-constraint (4C) energy-momentum conservation 
kinematic fit is performed to the J/tp K + K~tt + 7t~jj 
hypothesis, and only events with xic ^ ^00 are accepted. 
For events with more than two photon candidates, the 
combination with the minimum xlc I s selected. After 
the 4C fits, the ir° and 77 signal windows on the 77 
invariant mass distributions are defined in the ranges 
0.115 < m 77 < 0.150 GeV/c 2 and 0.518 < m 77 < 0.578 
GeV/c 2 , respectively. 

The numbers of J/tp — > <prj (77'), 77 (77') — > 
tt + ir~ ir° (77) events are obtained from an unbinncd 
extended maximum likelihood (ML) fit to the K + K~ 
versus ir + ir~ ir° (77) invariant mass distributions. The 
projection of the fit on the ttikk (jn^^^-^o and m^+^-jj) 
axis is shown in Figs. [3] (a) and|4] (a) [Figs. [3] (b) and [4] 
(b)] for the 77 and 77' cases, respectively. In the ML fits, 
we require that 0.99 GeV/c 2 < m K K < 1-09 GeV/c 2 and 
0.50 GeV/c 2 < m,!,-,. < 0.60 GeV/c 2 (0.87 GeV/c 2 
< m^+ 7V - r) < 1.07 GeV/c 2 ) for the 77 (77') case. The signal 
shape for <p is modeled with a rclativistic Breit-Wigncr 
(RBW) function [34| convoluted with a Gaussian 
function that represents the detector resolution. In 
the fit, the width of <j> is fixed at the PDG value, 
and its central mass value is floated, and the width 
of the Gaussian is free; the signal shape for 77 (77') is 
described by a Crystal Ball (CB) function [35j], and its 
parameters are floated. The backgrounds are divided into 
three categories: non-</>77(77')-pcaking background (i.e., 
J/tp ^ Tr + ir~ir°K + K~); non-</>peaking background 
[i.e., J/tp —± K + K~r] (77')]; and non-77(77')-pcaking 
background (i.e., J/tp — > (jm + ir~ir a ). The probability 
density functions (PDF) for non-0-peaking background 



in the tukk distribution is parameterized by [36j 

B(m KK ) = (m KK - 2m K ) a ■ e r hm ^ m (1) 

where a, b and c are free parameters, and mx is 
the nominal mass value of the charged kaon from 
the PDG [2^]. The shape for the non-77 ( 77')-peaking 
background in the m^+^-^oi^ distribution is modeled 
by a first-order Chebychev polynomial function 
[B(mK+K-K°(r)))\- All parameters related to the 
background shape are floated in the fit to data. Totally, 
14 parameters including signal and background yields are 
floated in the fit. The PDFs for signal and backgrounds 
are combined in the likelihood function C, defined as a 



function of the free parameters TV 1 ', TV, 



and Kkg"" 1 



non-07/ 
bkg ' 



TV, 



non-0 
bkg ' 



TV! 



x J[[N*RBW(rrj KK ) x CB^,.^) 
i=i 

+N™^B(m i KK )xB(mi +7r _ n0 ) 
+N™-*B(m i KK ) x CB(mi +n ^ ) 
+N™«-«RBW(m i KK ) x B^)], (2) 

where TV'' is the number of J/ip — > cprj, <p — > K + K~ , 77 — > 
^"Tr events, and N^, and N^> 

are the numbers of the corresponding three kinds of 
backgrounds. The fixed parameter TV is the total number 
of selected events in the fit region, and rn l KK (m l + _ ) 
is the value of ttlkk (m^+^-^o) for the ith event. We 
use the product of the PDFs, since we have verified that 
ttlkk and m^+^- n o are uncorrelated for each component. 
The negative log-likelihood (— ln£) is then minimized 
with respect to the extracted yields. The resulting 
fitted signal and background yields are summarized in 
Table [H We also obtain the results for the 77' case by 
replacing 77 and ir° with 77' and 77 in Eq. ([2]). The fitted 
results for 77 (77') — > 7r + ir~ ir° (77) are shown in Fig. [3] 
(Fig. 2J. The detection efficiencies are determined with 
MC simulations to be 20.37% and 20.89% for 77 and 77', 
respectively. 



IV. SYSTEMATIC UNCERTAINTIES 

Contributions to the systematic error on the ratios 
are summarized in Table [TXl The uncertainty, due to the 
requirement of no neutral showers in the EMC inside a 
cone spanning 0.3 (1.5) rad around the direction of the 
missing momentum for 77 (77') decay is obtained using 
the control sample of decays J/tp — > (prf ', <p — > K + K~, 
77' — > 77r + 7r~. We calculated the missing momentum 
of the K + K~tt + it~ system, and define the same cones 
around the direction of the missing momentum as in the 77 
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FIG. 3: The (a) itlkk and (b) m 7T + n - 7r o distributions with fit results superimposed for J/ip — > <j>r), <f> — » K + K~ , r\ — > tt + tt~it . 
Points with error bars are data. The (black) solid curves show the results of the total fits, and the (black) short-dashed 
curves are for signals. The (blue) dotted-dash curve shows non-77-peaking backgrounds, the (red) dashed curve shows the 
non-<^>-peaking background, and the (green) dotted-short-dash curve shows non-0?7-peaking backgrounds. 




FIG. 4: The (a) thkk and (b) m n +^- rl distributions with fit results superimposed for J/ip — ► <jyq' , <p — > K + K~ , rj' — > ir + iv~r]. 
Points with error bars are data. The (black) solid curves show the results of the total fits, and the (black) short-dashed 
curves are for signals. The (blue) dotted-dash curve shows non-r/-peaking backgrounds, the (red) dashed curve shows the 
non-0-peaking background, and the (green) dotted-short-dash curve shows non-<^r/-peaking backgrounds 



(7/) semilcptonic analysis. The ratios of events with the 
requirement on the number of extra photons to events 
without the requirement are obtained for both data and 
MC simulation. The difference 0.1% (1.1%) is considered 
as a systematic error for the 77 (rj') case. 

We also use the control sample of J/ip — > (fn)', <p —> 
K + K~ , rf — > 77r + 7T~ to obtain the uncertainty due to 
the requirement on the missing momentum P m i ss > 0.03 
GeV/c for both r\ and 77' cases. Thus we calculated the 
missing momentum of the K + K^tt + tt^ system. The 
ratios of events with the requirement on the missing 
momentum P m i ss > 0.03 GeV/c to events without 
the requirement are obtained for both data and MC 
simulation. The difference 2.5% is considered 
systematic error for both r\ and 77' cases. 

The phase space MC is used to generate 77 (if) —> 



ir + e~D e decays. In Rcf. [12J, the transition form factors 
f±* are calculated at the one-loop level in the chiral 
perturbation theory. We use the model predictions to 
generate signal MC events, and find that the uncertainty 
on the detection efficiency is changed by 1.0% (5.0%) for 
the 77 (77') case. 

Since the uncertainties on the PID of the electron 
and one of the pions do not cancel in the ratio, the 
efficiencies for pion and electron PID are obtained with 
the control samples of J/i/j —> ir + ir~ir a and radiative 
Bhabha scattering e + e~ — > 7e + e~ (including J/ip —> 
7e + e _ ), respectively. S amp les with backgrounds less 
than 1.0% arc obtained [33] - The differences between 
data and MC for the efficiencies of pion and electron PID 
are about 1.0% and 1.2%, respectively, which arc taken 
as systematic errors. Using the same control samples, 
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TABLE I: The fitted signal and background yields for J/ip — > 
4>V V (v') 7r + 7r _ 7r°(ry), and € rl (e v ) is its selection 

efficiency. 



Quantity 


Value 
77 77' 


N V (N V ') 


3850 ± 73 1623 ±44 


» r non-</>77/ »rHOIl-^t)' \ 
Jv bkg V iv bkg ) 


24 ± 8 49 ± 10 


» riion-^ / iuiion-^i 
Jv bkg ^ v bkg J 


367 ±43 22 ± 17 


jy-non-Tj/iy-non-V', 
Jv bkg l Jv bkg J 


88 ±14 61 ± 12 


£"(6"') 


20.37% 20.89% 



we estimate the uncertainty due to the requirement of 
E/p for the electron selection to be 3.5% (3.4%) for the 
77 (77') case, and the uncertainty due to the requirement 
of E/p for pion selection in the 77 semileptonic decay is 
estimated to be 0.8%. The systematic uncertainty due 
to the requirements of <f> and 77 (77') mass windows are 
estimated to be 1.4% and 0.04% (0.2%) by using the 
control sample of J/ijj — > (fnj (77'), 77 (77') — > ir + tt~ tt° (77) . 



The uncertainty in the determination of the numbers 
of observed events for J/ip — > (prj [77 — > 7r + 7r~7r 
77)] and J/ip — » 4>f]' [77' — > ir + n~ri (77 — > 77)] 
decays are estimated on the basis of earlier published 
results. The photon detection efficiency and its 
uncertainty are studied by three different methods in 
Ref [33 ■ The systematic error of photon detection is 
estimated to be 1.0% per photon. In the fit to the 
</> mass distribution, the mass resolution is fixed to 
the MC simulation; the level of possible discrepancy 
is determined with a smearing Gaussian, for which a 
non-zero a would represent a MC-data difference in 
the mass resolution. The uncertainty associated with 
the difference determined in this way is 0.03% (0.06%) 
for the 77 (77') case. The systematic uncertainty due 
to the choice of parameterization for the shape of the 
non-</>77 ( 77')-peaking background is estimated by varying 
the order of the polynomial in the fit; we find the relative 
changes on the 77 (77') signal yield of 1.3% (0.8%), which is 
taken as the uncertainty due to the background shapes. 
The systematic errors from ir° (77) reconstruction from 
77 decays is determined to be 1.0% per ir° (77) by using a 
high purity control sample of J/ip — > ir°pp (J/ip ijpp) 
decay [HI. The branching fractions for the tt° and 
77 — >• 77 decays arc taken from the PDG [29|. The 
uncertainties on the branching fractions are taken as a 
systematic uncertainty in our measurements. The total 
systematic error a^ ys (cr^/ s ) on the ratio is 5.6% (7.4%) 
for 77 (77'), as summarized in Table ITI1 



TABLE II: Summary of relative systematic errors for the 
determination of ratios of semileptonic over hadronic decays. 
The first nine lines are relevant for the semileptonic weak 
decay chain J/ip — > (prj (77'), 77 (77') — ► ir e~v e . The next 
five lines are relevant for the determination of the signal yield 
of the hadronic decay process J/ip — > cpV V ~^ t{ + tv~-it 
(t/ — > ir + ir~r]). 





Sys. 


error (%) 


Sources 


V 




Requirement on iVshower 


0.1 


1.1 


Requirement on P m i ss 


2.5 


2.5 


Signal model 


1.0 


5.0 


Electron PID 


1.2 


1.2 


Requirement on E/p for e 


3.5 


3.4 


Requirement on E/p for n 


0.8 




4> mass window 


1.4 


1.4 


77 (77') mass window 


0.0 


0.2 


Photon efficiency 


2.0 


2.0 


7T PID 


1.0 


1.0 


Signal shapes for 77(77') — > 7r + 7r _ 7r (77) 


0.0 


0.1 


Background shape for 77(77') — s> 7r + 7r~7r (7)) 


1.3 


0.8 


7r°(?7) reconstruction 


1.0 


1.0 


Cited /3(tt (77) -^77) 


0.0 


0.5 


Total 


5.6 


7.4 



V. RESULTS 

The upper limit on the ratio of branching fractions of 
the semileptonic decay B(rj ir + e~i> e + c.c.) over the 
hadronic decay B(r\ — > tt + tt~tt°) is calculated with 



73(77 — > ir + e v e ± c.c.) 
B(rj — > ir + 7T~ir°) 



< 



77) 



Nv/e* (1 - <t„) 



(3) 



where N^ L is the 90% upper limit of the observed number 
of events for J/ip — > <pr/ 7 <p —> K + K~, 77 — > Tr + e~v e 
decay, e^ L is the MC determined efficiency for the signal 
channel, N n is the number of events for the J/ip — > cpi], 
<p — > K + K~, 77 -4 7r + 7r~7r°(7r° — » 77) decay, e v is 
the MC determined efficiency for the decay mode, and 



a 



5.9%, where a* ys and cr" tat 



stat ^2 

77 — y \"V 1 T V u T f ) 

are the total relative systematic error for the 77 case 
from Table [II] and the relative statistical error of N v , 
1 



7.9% 



0. 



respectively. For 77', cy = J{<r a j s ) 2 + (<^ at ) 2 

The relative statistical error of TV 1 ' (iV') is 1.9% (2.7%). 
We also obtain the upper limit on the ratio of B(rj' — > 
TT + e~D e + c.c.) to B(rf — > 7r + 7r~?7) by replacing 77 and 
B(n° — > 77) with 77' and B(r\ —¥ 77), respectively, in 
Eq. ((3]). Since only the statistical error is considered 
when we obtain the 90% upper limit of the number of 

events, to be conservative, N^ L and N^ L are shifted up 
by one sigma of the additional uncertainties (<j n or ay). 
Using the numbers in Table IIII1 the upper limits on 

tViP rating g (')-> 7r + e ~ Pe+c.c.) , B(rt' e~ v c +c.c. 

trie ratios 8 ( v ^+ir-^) and — B( n '^ir+ir- v ) 
obtained at the 90% C.L. of 7.3 x 10~ 4 and 5.0 x 10~ 4 
respectively. 



arc 
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TABLE III: The numbers used in the calculations of the ratios 
in Eq. ©, where ^Vul(-^ul) ^ s * ne u PP er limit of the signal 
events at the 90% C.L., e^ L (e^) is the selection efficiency, 

N V (N V ) is the number of the events of J/ip — > <f>r]{ri'), <f) — > 
K + K~ , r\ — > TT + ir~TT° and ir° —¥ 77 (77' — > TT + ir~ri and r\ —¥ 
77) , e' 7 (e v> ) is its selection efficiency, cr® tat (a^ at ) is the 

relative statistical error of N v (N v ) , and a v (<t, ; < ) is the total 
relative error. 



Quantity 


Value 
r? 77' 


N v (N v ') 
e" (e"') 

(7,, (<7y) 


2.36 1.59 
17.9% 17.4% 
3850 ± 73 1623 ±44 
20.37% 20.89% 
1.9% 2.7% 
5.9% 7.9% 



VI. SUMMARY 

A search for the semileptonic weak deacy modes 
V iv') ~* 7i" + e _ P e has been performed for the first time in 
the process of J/ip ->■ 077(77') using the (225.3 ±2.8) x 10 6 
J/ip events measured at BESIII. We find no signal 
yields for the semileptonic weak decays of 77 and 77'. 
The upper limits at the 90% C.L. arc 7.3 x 1CT 4 and 
5.0 x 10~ 4 for the ratios of semileptonic over hadronic 

decay modes s^^+^o) and — B ( v '^+-n- v ) — > 
respectively. The advantage of measuring the ratios 
instead of the branching fractions of the semileptonic 
weak decays is that many uncertainties cancel. Using 
the hadronic branching fraction values of 77 — > 7r+7r~7r° 



and 77' — ► 7T+7T 77 as listed by PDG (29[, we determine the 
semileptonic decay rates to be B(r\ — > 7r + e~P e + c.c.) < 
1.7 x 10~ 4 and 8(77' -> 7r+e^ e ± c.c.) < 2.2 x 10~ 4 at 
the 90% C.L.. 
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